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Abstract. A pedagogical derivation is presented of the "fireball" model of 7-ray 
bursts, according to which the observable effects are due to the dissipation of the ki- 
netic energy of a relativistically expanding wind, a "fireball." The main open questions 
are emphasized, and key afterglow observations, that provide support for this model, 
are briefly discussed. The relativistic outflow is, most likely, driven by the accretion of 
a fraction of a solar mass onto a newly born (few) solar mass black hole. The observed 
radiation is produced once the plasma has expanded to a scale much larger than that 
of the underlying "engine," and is therefore largely independent of the details of the 
progenitor, whose gravitational collapse leads to fireball formation. Several progenitor 
scenarios, and the prospects for discrimination among them using future observations, 
are discussed. The production in 7-ray burst fireballs of high energy protons and neu- 
trinos, and the implications of burst neutrino detection by kilometer-scale telescopes 
under construction, are briefly discussed. 



1 Introduction 

The widely accepted interpretation of the phenomenology of 7-ray bursts (GRBs), 
of 0.1 — 1 MeV photons lasting for a few seconds (see for a review), is that the 
observable effects are due to the dissipation of the kinetic energy of a relativis- 
tically expanding wind, a "fireball," whose primal cause is not yet known. The 
recent detection of "afterglows," delayed low energy (X-ray to radio) emission 
of GRBs (see jSHj for a review) , confirmed the cosmological origin of the bursts 
through the rcdshift determination of several GRB host-galaxies, and supported 
standard model predictions of afterglows that result from the collision of an 
expanding fireball with its surrounding medium (see |78I92) for reviews). 

The fireball model is described in §§1201 and 0] of this chapter. The phe- 
nomcnological arguments suggesting that fireball formation is likely regardless of 
the underlying progenitor are presented, and fireball hydrodynamics and radia- 
tive processes are discussed in detail in §33 and [21 respectively. The main open 
questions related to fireball physics are discussed in M3.4I Since both the theory 
and the implications of afterglow observations are extensively discussed in other 
chapters of this volume, we include in 21 of this chapter only a brief discussion 
of several key afterglow implications. We also limit the theoretical discussion of 
fireball evolution to the GRB production phase that precedes the afterglow phase 
during which evolution is dominated by the interaction of the fireball with its 
surrounding medium. We do discuss, however, the initial non-self-similar onset 
of this interaction, which marks the onset of the afterglow phase. 
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The GRB progenitors are not yet known. We present in SjSJthe constraints 
imposed by observations on possible progenitors, and discuss the (presently) 
leading candidates. Hints provided by afterglow observations, which are exten- 
sively discussed in separate chapters of this volume, are briefly reviewed. 

The association of GRBs with ultra- high energy cosmic-rays (UHECR), the 
evidence for which is strengthened by recent afterglow observations, is based on 
two key points |115| : 1) the constraints that a dissipative ultra-relativistic wind 
must satisfy in order to allow acceleration of protons to energy ~ 10 20 eV, the 
highest observed cosmic-ray energy, are remarkably similar to the constraints 
imposed on the fireball wind by 7-ray observations, and 2) the inferred local 
(z = 0) GRB energy generation rate of 7-rays is remarkably similar to the local 
generation rate of UHECR implied by cosmic-ray observations. We briefly discuss 
in [j^l production of high energy protons and neutrinos in GRB fireballs (see 
[12()I121| for more detailed reviews). The GRB model for UHECR production 
makes unique predictions that may be tested with operating and planned large 
area UHECR detectors |27l2yil()filll4| 1 . In this review we focus, however, on 
more recent predictions of neutrino emission, which may be tested with planned 
high energy neutrino telescopes |S7]. Detection of the predicted neutrino signal 
will confirm the GRB fireball model for the production of UHECR and may 
allow discrimination between different fireball progenitor scenarios. Moreover, a 
detection of even a handful of neutrino events correlated with GRBs will allow a 
test for neutrino properties, e.g., flavor oscillation and coupling to gravity, with 
accuracy many orders of magnitude better than currently possible. 

2 The Fireball Model: Hydrodynamics 
2.1 Relativistic Expansion 

General phenomenological considerations, based on 7-ray observations, indicate 
that, regardless of the nature of the underlying sources, GRBs are produced 
by the dissipation of the kinetic energy of a relativistically expanding fireball. 
The rapid rise time and short duration, ~ 1 ms, observed in some bursts [T3 38 
imply that the sources are compact, with a linear scale comparable to a light- 
ms, To ~ 10 cm. The high 7-ray luminosity implied by cosmological distances, 
L 1 ~ 10 52 erg s _1 , then results in a very high optical depth to pair creation since 
the energy of observed 7-ray photons is above the threshold for pair production. 
The number density of photons at the source n 7 is approximately given by 

L-y = 47rrgcn 7 e, (1) 

where e ~ lMeV is the characteristic photon energy. Using tq ~ 10 7 cm, the 
optical depth for pair production at the source is 

1 See also http://www.physicsMdelaide.edu.au/astrophysics/FlysEye.html, 
http://www-ta.icrr.u-tokyo.ac.jp/, and http://www.auger.org/ 
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where ctt is the Thomson cross section. 

The high optical depth implies that a thermal plasma of photons, electrons, 
and positrons is created, a fireball which then expands and accelerates to rela- 
tivists velocities [50187) . The optical depth is reduced by relativistic expansion of 
the source. If the source expands with a Lorentz factor r, the energy of photons 
in the source frame is smaller by a factor r compared to that in the observer 
frame, and most photons may therefore be below the pair production threshold. 

A lower limit for F may be obtained in the following way [10167) . The GRB 
photon spectrum is well fitted in the Burst and Transient Source Experiment 
(BATSE) detectors range, 20 keV to 2 MeV [35], by a combination of two power- 
laws, dn~/de„ oc e 7 " (a is the flux density spectral index, F v cx v +a ) with 
different values of a at low and high energy [5]. Here, dn 1 /de 1 is the number 
of photons per unit photon energy. The break energy (where a changes) in the 
observer frame is typically e 7 b ~ 1 MeV, with a ~ at energies below the break 
and a ~ — 1 above the break. In several cases, the spectrum has been observed to 
extend to energies > 100 MeV Consider then a high energy test photon, 

with observed energy e t , trying to escape the relativistically expanding source. 
Assuming that, in the source rest frame, the photon distribution is isotropic, 
and that the spectrum of high energy photons follows dn~ 1 /de 1 oc e~ 2 , the mean 
free path for pair production (in the source rest frame) for a photon of energy 
e' t = et/7 (in the source rest frame) is 



Here, e t h {e't , 0) is the minimum energy of photons that may produce pairs inter- 
acting with the test photon, given by e^e^l — cos 6*) > 2(m e c 2 ) 2 (9 is the angle 
between the photons' momentum vectors). U 1 is the photon energy density (in 
the range corresponding to the observed BATSE range) in the source rest-frame, 
given by 



(Note that we have used a constant cross section, 3<7t/16, above the threshold 
eth-) The cross section drops as log(e)/e for e 3> eth; however, since the number 
density of photons drops rapidly with energy, this does not introduce a large 
correction to i 77 . 

The source size constraint implied by the variability time is modified for a 
relativistically expanding source. Since in the observer frame almost all photons 
propagate at a direction making an angle < l/T with respect to the expansion 
direction, radiation seen by a distant observer originates from a conical section 
of the source around the source-observer line of sight, with opening angle ~ 1/-T. 
Photons which are emitted from the edge of the cone are delayed, compared to 
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those emitted on the line of sight, by r/2r 2 c. Thus, the constraint on source size 
implied by variability on time scale At is 



2r 2 cAt. 



(5) 



The time r/c required for significant source expansion corresponds to comoving 
time (measured in the source frame) i co . w r/Tc. The two-photon collision rate at 

c/ljy. Thus, the source optical depth to pair production 
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n /ri 71 . Using Eqs. © and JHJ) we have 
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r > 250 



L 



7,52 



( ^ ) 

VIOOMeV/ 



where L~ 



10 52 i 7 , 52 erg s _1 and At 



lOOMeV 
\Q- 2 At 



Atzl 



1/6 



(6) 



(7) 



-2 S. 



2.2 Fireball Evolution 

As the fireball expands it cools, the photon temperature T 7 in the fireball frame 
decreases, and most pairs annihilate. Once the pair density is sufficiently low, 
photons may escape. However, if the observed radiation is due to photons escap- 
ing the fireball as it becomes optically thin, two problems arise. First, the photon 
spectrum is quasi-thermal, in conflict with observations. Second, the source size, 
ro ~ 10 7 cm, and the total energy emitted in 7-rays, ~ 10 53 erg, suggest that the 
underlying energy source is related to the gravitational collapse of a ~ 1 M 
object. Thus, the plasma is expected to be loaded with baryons which may be 
injected with the radiation or present in the atmosphere surrounding the source. 
A small baryonic load, > 10 -8 M , increases the optical depth due to Thomson 
scattering on electrons associated with the loading protons, so that most of the 
radiation energy is converted to kinetic energy of the relativistically expanding 
baryons before the plasma becomes optically thin |88ll()2j . To overcome both 
problems it was proposed jHS] that the observed burst is produced once the ki- 
netic energy of the ultra-relativistic ejecta is re-randomized by some dissipation 
process at large radius, beyond the Thomson photosphere, and then radiated as 
7-rays. Collision of the relativistic baryons with the interstellar medium (ISM) 
|95| . and internal collisions within the ejecta itself |79I86I91] . were proposed as 
possible dissipation processes. Most GRBs show variability on time scales much 
shorter (typically 10 -2 times) than the total GRB duration. Such variability 
is hard to explain in models where the energy dissipation is due to external 
shocks 98 127 . Thus, it is believed that internal collisions are responsible for 
the emission of 7-rays. 

Let us first consider the case where the energy release from the source is 
"instantaneous," i.e., on a time scale of rvj/c. We assume that most of the energy 
is released in the form of photons, i.e., that the fraction of energy carried by 
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baryon rest mass M satisfies rj^ 1 = Mc 2 /E <C 1. The initial thickness of the 
fireball shell is tq. Since the plasma accelerates to relativistic velocity, all fluid 
elements move with velocity close to c, and the shell thickness remains constant 
at r (this changes at very late time, as discussed below). We are interested in 
the stage where the optical depth (due to pairs and/or electrons associated with 
baryons) is high, but only a small fraction of the energy is carried by pairs. 

The entropy of a fluid component with zero chemical potential is S — V (e + 
p)/T, where e, p and V are the (rest frame) energy density, pressure and vol- 
ume. For the photons p = e/3 oc Since initially both the rest mass and 
thermal energy of baryons is negligible, the entropy is provided by the photons. 
Conservation of entropy implies 

r 2 r(r)r a T^(r) = const, (8) 
and conservation of energy implies 

r 2 r(r)r r(r)T*(r) = const . (9) 
Here r(r) is the shell Lorentz factor. Combining (JSJ) and 10 we find 

r(r) oc r, Tj oc r^ 1 , n oc r -3 , (10) 

where n is the rest frame (comoving) baryon number density. 

As the shell accelerates the baryon kinetic energy, rMc 2 , increases. It be- 
comes comparable to the total fireball energy when r ~ rj, at radius rf = rjro- 
At this radius most of the energy of the fireball is carried by the baryon kinetic 
energy, and the shell does not accelerate further. Eq. JHJ describing energy con- 
servation is replaced with r = const. Eq. ijHJl, however, still holds. Eq. |(HJ may be 
written as / nT 7 = const (constant entropy per baryon) . This implies that the 
ratio of radiation energy density to thermal energy density associated with the 
baryons is r independent. Thus, the thermal energy associated with the baryons 
may be neglected at all times, and Eq. JHJ holds also for the stage where most 
of the fireball energy is carried by the baryon kinetic energy Thus, for r > rf 
we have 

r(r)=r&r), Txr- 2/3 , uocr 2 . (11) 

Let us consider now the case of extended emission from the source, on time 
scale 3> tq/c. In this case, the source continuously emits energy at a rate L, and 
the energy emission is accompanied by mass-loss rate 

M = L/?]c 2 . (12) 

For r < rf the fluid energy density is relativistic, aT^/nm p c 2 = r/r^/r, and 
the speed of sound is ~ c. The time it takes the shell at radius r to expand 
significantly is r/c in the observer frame, corresponding to t co . ~ r/Tc in the 
shell frame. During this time sound waves can travel a distance cr/Tc in the shell 
frame, corresponding to r / T 2 — r / (r/rp) 2 = (ro/ r ) r o m the observer frame. This 
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implies that at the early stages of evolution, r ~ r , sound waves had enough 
time to smooth out spatial fluctuations in the fireball over a scale ro, but that 
regions separated by Ar > ro cannot interact with each other. Thus, if the 
emission extends over a time ?;grb 3> r /c, a fireball of thickness c£qrb S> r 
would be formed, which would expand as a collection of independent, roughly 
uniform, sub-shells of thickness r . Each sub-shell would reach a final Lorentz 
factor i~y, which may vary between sub-shells. This implies that different sub- 
shells may have velocities differing by Av ~ c/2?7 2 , where 77 is some typical value 
representative of the entire fireball. Different shells emitted at times differing 
by At, r n /c < At < £grb, may therefore collide with each other after a time 
t c ~ cAt/Av, i.e., at a radius 

n w 2r 2 cAt = 6 x 10 13 r 2 2 5 Z\t_ 2 cm, (13) 

where r = lO 2 5 /^. The minimum internal shock radius, r ~ -T 2 r , is also the 
radius at which an individual sub-shell may experience significant change in its 
width ro, due to Lorentz factor variation across the shell. 

2.3 The Allowed Range of Lorentz Factors and Baryon Loading 

The acceleration, r oc r, of fireball plasma is driven by radiation pressure. Fire- 
ball protons are accelerated through their coupling to the electrons, which are 
coupled to fireball photons. We have assumed in the analysis presented above, 
that photons and electrons are coupled throughout the acceleration phase. How- 
ever, if the baryon loading is too low, radiation may decouple from fireball elec- 
trons already at r < rj. The fireball Thomson optical depth is given by the 
product of comoving expansion time, r/r(r)c, and the photon Thomson scatter- 
ing rate, n e ccr<r. The electron and proton comoving number densities arc equal, 
n e = rip, and are determined by equating the r independent mass flux carried 
by the wind, 4:Trr 2 cr(r)n p m p , to the mass-loss rate from the underlying source, 
which is related to the rate L at which energy is emitted through M = L/(ric 2 ). 
Thus, during the acceleration phase, r(r) = r/r n and the Thomson optical 
depth, tt, is tt oc r~ 3 . The Thomson optical depth drops below unity at a 
radius r < r/ = -qr^ if r\ > 77*, where 

Here ro = 10 7 ro,7 cm. 

If 77 > 77* radiation decouples from the fireball plasma at r = r/ro = 
ijt If r) 77*, then most of the radiation energy is not converted to 

kinetic energy prior to radiation decoupling, and most of the fireball energy 
escapes in the form of thermal radiation. Thus, the baryon load of fireball 
shells, and the corresponding final Lorentz factors, must be within the range 
10 2 < r w 77 < 77* w 10 3 in order to allow the production of the observed 
non-thermal 7-ray spectrum. 
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2.4 Fireball Interaction with the Surrounding Medium 

As the fireball expands, it drives a relativistic shock (blastwave) into the sur- 
rounding gas, e.g., into the ISM gas if the explosion occurs within a galaxy. In 
what follows, we refer to the surrounding gas as ISM gas, although the gas need 
not necessarily be interstellar. At early times, the fireball is little affected by 
the interaction with the ISM. At late times, most of the fireball energy is trans- 
ferred to the ISM, and the flow approaches the self-similar blastwave solution 
of Blandford and McKee [T2|. At this stage a single shock propagates into the 
ISM, behind which the gas expands with Lorentz factor 

r B u(r) = ( ~rF~~~ 2 ^ V2 ^ 3/2 = 150 r~r\ (15) 

\lbTTnm p c z J \ n a J 

where E = 10 53 £ , 53 erg is the fireball energy, n = Iuq cm -3 is the ISM number 
density, and r = 10 17 rn cm is the shell radius. The characteristic time at which 
radiation emitted by shocked plasma at radius r is observed by a distant observer 
is t ks r/4r BM c [IIS]. 

The transition to self-similar expansion occurs on a time scale tss (measured 
in the observer frame) comparable to the longer of the two time scales set by the 
initial conditions: the (observer) GRB duration £grb and the (observer) time tr 
at which the self-similar Lorentz factor equals the original ejecta Lorentz factor 
r, r BM (t = t r ) = r. Since t = r/Arl M c, 



t = max 



, r , ^53 \ ^ 3 -8/3 

iGRB,5 ( — J r 25 > s 



(16) 



During the transition, plasma shocked by the reverse shocks expands with Lorentz 
factor close to that given by the self-similar solution, 

/ rp \i/8 

r tr . ~ r BM (t = i ss ) ~ 245 f J ^ 3/8 , (17) 

where t — 10ii s. The unshocked fireball ejecta propagate at the original expan- 
sion Lorentz factor, r, and the Lorentz factor of plasma shocked by the reverse 
shock in the rest frame of the unshocked ejecta is ~ r/r tv .. If t ~ £grb ~S> tr 
then r/r tI . ^> 1, the reverse shock is relativistic, and the Lorentz factor associ- 
ated with the random motion of protons in the reverse shock is ~ f/ftr.. 

If, on the other hand, t ~ tr 3> ^grb then r/Tu. ~ 1, and the reverse 
shock is not relativistic. Nevertheless, the following argument suggests that the 
reverse shock speed is not far below c, and that the protons are therefore heated 
to relativistic energy, — 1 ~ 1. The comoving time, measured in the fireball 
ejecta frame prior to deceleration, is t co ~ r/Tc. The expansion Lorentz factor is 
expected to vary across the ejecta, Ar/T ~ 1, due to variability of the underlying 
GRB source over the duration of its energy release. Such variation would lead 
to expansion of the ejecta, in the comoving frame, at relativistic speed. Thus, 
at the deceleration radius, t co . ~ rt, the ejecta width exceeds ~ rf co . ~ Fct. 
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Since the reverse shock should cross the ejecta over a deceleration time scale, 
~ rt, the reverse shock speed must be close to c. We therefore conclude that the 
Lorentz factor associated with the random motion of protons in the reverse shock 
is approximately given by — 1 ~ r/r tr . for both r/r tr . ~ 1 and r/r tI . 3> 1. 

Since icRB ^ 10 s is typically comparable to tp, the reverse shocks are 
typically expected to be mildly relativistic. 

2.5 Fireball Geometry 

We have assumed in the discussion so far that the fireball is spherically sym- 
metric. However, a jet-like fireball behaves as if it were a conical section of a 
spherical fireball as long as the jet opening angle is larger than _T _1 . This is due 
to the fact that the linear size of causally connected regions, ct co . ~ r/T in the 
fireball frame, corresponds to an angular size ct co ./r ~ r . Moreover, due to 
the relativistic beaming of radiation, a distant observer cannot distinguish be- 
tween a spherical fireball and a jet-like fireball, as long as the jet opening angle 
6 > r . Thus, as long as we are discussing processes that occur when the wind 
is ultra-relativistic, r ~ 300 (prior to significant fireball deceleration by the sur- 
rounding medium), our results apply for both a spherical and a jet-like fireball. 
In the latter case, L (E) in our equations should be understood as the luminosity 
(energy) the fireball would have carried had it been spherically symmetric. 

3 The Fireball Model: Radiative Processes 
3.1 Gamma- Ray Emission 

If the Lorentz factor variability within the wind is significant, internal shocks will 
reconvert a substantial part of the kinetic energy to internal energy. The internal 
energy may then be radiated as 7-rays by synchrotron and inverse Compton 
emission of shock-accelerated electrons. The internal shocks are expected to be 
mildly relativistic in the fireball rest frame, i.e., characterized by Lorentz factor 
.Tj — 1 ~ a few. This is due to the fact that the allowed range of shell Lorentz 
factors is ~ 10 2 to ~ 10 3 (see ^2.3(1 . implying that the Lorentz factors associated 
with the relative velocities are not very large. Since internal shocks are mildly 
relativistic, we expect results related to particle acceleration in sub-relativistic 
shocks (see ^B] for a review) to be valid for acceleration in internal shocks. 
In particular, electrons are expected to be accelerated to a power law energy 
distribution, dn e /dr e oc r~ p for T e > r m , with p ~ 2 |5I12I18| . 

The minimum Lorentz factor r m is determined by the following considera- 
tion. Protons are heated in internal shocks to random velocities (in the wind 
frame) I^f — 1 wTj- If electrons carry a fraction £ e of the shock inter- 

nal energy, then F m « £ e (m p /m e ). The characteristic frequency of synchrotron 
emission is determined by F m and by the strength of the magnetic field. As- 
suming that a fraction £b of the internal energy is carried by the magnetic field, 
Anrf cr 2 B 2 / 8tt = £s-Li n t., the characteristic observed energy of synchrotron pho- 
tons, e 7 fe = rhr^eB/rrieC, is 
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e jb » MeV. (18) 

In deriving Eq. I|18l) we have assumed that the wind luminosity carried by internal 
plasma energy, Lint.) is related to the observed 7-ray luminosity through L mt , — 
L 7 /£ e . This assumption is justified since the electron synchrotron cooling time 
is short compared to the wind expansion time (unless the equipartition fraction 
£b is many orders of magnitude smaller than unity), and hence electrons lose 
all their energy radiatively. Fast electron cooling also results in a synchrotron 
spectrum c?n 7 /c?e 7 oc e 7 1 p ^ 2 — e~ 2 at e 7 > e 7 h, consistent with observed GRB 
spectra [§]. 

At present, there is no theory that allows the determination of the values 
of the equipartition fractions £ e and £b. Eq. (|18l) implies that fractions not far 
below unity are required to account for the observed 7-ray emission. We note that 
build up of magnetic field to near equipartition by electromagnetic instabilities is 
expected to be a generic characteristic of collisionless shocks (see the discussion 
in |16| and references therein) , and is inferred to occur in other systems such as 
in supernova remnant shocks (see, e.g., |24l6(Jp . 

3.2 Break Energy Distribution 

The 7-ray break energy e lb of most GRBs observed by BATSE detectors is in 
the range of 100 keV to 300 keV [22]. It may appear from Eq. (|T%j> that the 
clustering of break energies in this narrow energy range requires fine tuning of 
fireball model parameters, which should naturally produce a much wider range of 
break energies. This is, however, not the case [55]. Consider the dependence of e 7 b 
on r. The strong r dependence of the pair production optical depth, Eq. (0, 
implies that if the value of r is smaller than the minimum value allowed by 
Eq. 10, for which T 77 (e 7 — lOOMeV) « 1, most of the high energy photons in the 
power-law distribution produced by synchrotron emission, <in 7 /de 7 cx e^ 2 , would 
be converted to pairs. This would lead to high optical depth due to Thomson 
scattering on e , and hence to strong suppression of the emitted flux |55) . For 
fireball parameters such that r 77 (e 7 = lOOMeV) ~ 1, the break energy implied 
by Eqs. JTBJ and Q is 

« l# a # 9 ^ MeV. (19) 
At 2 !* 

As explained in shell Lorentz factors cannot exceed 77* ~ 10 3 , for which 

break energies in the X-ray range, e lb ~ 10 keV, may be obtained. We note, 
however, that the radiative flux would be strongly suppressed in this case too 
|55| . If the typical F of radiation emitting shells is close to 77*, then the range 
of Lorentz factors of wind shells is narrow, which implies that only a small 
fraction of wind kinetic energy would be converted to internal energy which can 
be radiated from the fireball. 
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Thus, the clustering of break energies at ~ 1 MeV is naturally accounted 
for, provided that the variability time scale satisfies At < 10~ 2 s, which implies 
an upper limit on the source size, since At > tq/c (see jH6l49] for alternative 
explanations). We note, that a large fraction of bursts detected by BATSE show 
variability on the shortest resolved time scale, ~ 10 ms |127) . In addition, a 
natural consequence of the model is the existence of low luminosity bursts with 
low, 1 — 10 keV, break energies Such "X-ray bursts" may have recently 
been identified |64| . 

For internal collisions, the observed 7-ray variability time, ~ ri/r 2 c ~ At, 
reflects the variability time of the underlying source, and the GRB duration 
reflects the duration over which energy is emitted from the source. Since the 
wind Lorentz factor is expected to fluctuate on time scales ranging from the 
shortest variability time rvj/c to the wind duration ^grBj internal collisions will 
take place over a range of radii, r ~ r 2 rg to r ~ -T 2 ciGRB- 

3.3 Afterglow Emission 

Let us consider the radiation emitted from the reverse shocks during the tran- 
sition to self-similar expansion. The characteristic electron Lorentz factor (in 
the plasma rest frame) is T m ~ ^ e (r/r tl -.)m p /m e , where the internal energy 
per proton in the shocked ejecta is ~ (r/r tr .)m p c 2 . The energy density U is 
E w 4irr 2 ctr 2 r U, and the number of radiating electrons is 7V e w E/rm p c 2 . Us- 
ing Eq. (|17fl and r — 4_T t 2 . ct, the characteristic (or peak) energy of synchrotron 
photons (in the observer frame) is |124| 

e im « nr tr rl^- c = 2^_ 1 ^/ 2 _X /2 A 2 .5 eV, (20) 

and the specific luminosity, L e = dL/de 7 , at e 7m is 

L m « (2 7 ra)- 1 r t ,^AT e « lO^e^j^V 74 ^ 1 ^ 4 "- 1 , (21) 

where £ e = 0.1^ e _i, and £ B = 0.1^,-1- 

Here too, we expect a power law energy distribution, dN e /dr e cx r~ p for 
r e > r m , with p ~ 2. Since the radiative cooling time of electrons in the reverse 
shock is long compared to the ejecta expansion time, the specific luminosity 
extends in this case to energy e 7 > e 7TO as L t = £ m (e 7 /e 7m ) -1 / 2 , up to photon 
energy e 7C . Here e 7C is the characteristic synchrotron frequency of electrons for 
which the synchrotron cooling time, 6Trm e c/aTjeB 2 , is comparable to the ejecta 
(rest frame) expansion time, ~ r tl j. At energy e 7 > e 7C , 

e 7C w O.l^ln^E^ 2 ^ 1/2 keV, (22) 

the spectrum steepens to L e oc e~ . 

The shock driven into the ISM continuously heats new gas, and produces 
relativistic electrons that may produce the delayed afterglow radiation observed 
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on time scales t iss> typically of order days to months. As the shock wave 
decelerates, the emission shifts to lower frequency with time. Since afterglow 
emission on such long time scale is extensively discussed in other chapters of this 
volume, we do not discuss in detail the theory of late-time afterglow emission. 

3.4 Open Questions: Magnetic Field and Electron Coupling 

The emission of radiation in both the GRB and afterglow phases is assumed to 
arise from synchrotron emission of shock accelerated electrons. To match obser- 
vations, the magnetic field behind the shocks must be close to equipartition and 
a significant fraction of the internal shock energy must be carried by electrons, 
that is, £b and £ e should be close to unity, of order 10%. During the afterglow 
phase, shock compression of the existing ISM field yields a field many orders of 
magnitude smaller than needed. Thus, the magnetic field is most likely generated 
in, and by, the shock wave. A similar process is likely necessary to generate the 
field required for synchrotron emission during the GRB phase, i.e., in the inter- 
nal fireball shocks. Although a magnetic field close to equipartition at the base 
of the wind frozen into the fireball plasma may not be many orders of magnitude 
below equipartition during the internal shock phase, significant amplification is 
nevertheless required. It is well known that near equipartition magnetic fields 
may be generated in collisionless shocks through the Weibel instability (see, 
e.g., [ESEZj)- However, the field is generated on microscopic, skin-depth, scale 
and is therefore expected to rapidly decay, unless its coherence length grows to 
a macroscopic scale [53154) . The process by which such scale increase is achieved 
is not understood, and probably related to the process of particle acceleration 

In order to produce the observed spectrum during both afterglow and GRB 
phases, electrons must be accelerated in the collisionless shocks to a power-law 
distribution, dn e /dr e oc r~ p with p ~ 2. As mentioned in b|3 . II such distri- 
bution is expected in the internal shocks, which are mildly relativistic. Recent 
numeric and analytic calculations of particle acceleration via the first order Fermi 
mechanism in relativistic shocks show that similar spectral indices, p ~ 2.2, are 
obtained for highly-relativistic shocks as well [11165] . The derivation of electron 
spectral indices is based, in both the non-relativistic and relativistic cases, on a 
phenomenological description of electron scattering and, therefore, does not pro- 
vide a complete basic principle description of the process. In particular, these 
calculations do not allow one to determine the fraction of energy £ e carried by 
electrons. 

4 The Fireball Model: Key Afterglow Implications 

The following point should be clarified in the context of afterglow observations - 
the distribution of GRB durations is bimodal, with broad peaks at £grb ~ 0.2 s 
and £grb ~ 20 s 39,. The majority of bursts belong to the long duration, 
^grb ~ 20 s, class. The detection of afterglow emission was made possible thanks 
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to the accurate GRB positions provided on hour time scale by the BeppoSAX 
satellite [28]. Since the detectors on board this satellite trigger only on long 
bursts, afterglow observations are not available for the smaller population of 
short, £grb ~ 0.2 s, bursts. Thus, while the discussion of the fireball model pre- 
sented in §32Jand[31 based on 7-ray observations and on simple phcnomcnological 
arguments, applies to both long and short duration bursts, the discussion below 
of afterglow observations applies to long duration bursts only. It should, there- 
fore, be kept in mind that short duration bursts may constitute a different class 
of GRBs which, for example, may be produced by a different class of progenitors 
and may have a different redshift distribution than the long duration bursts. 

Afterglow observations led to the confirmation, as mentioned in §1, of the 
cosmological origin of GRBs and supported |117I126| standard model pre- 
dictions 62 80 90 112 of afterglows that result from synchrotron emission by 
electrons accelerated to high energy in the highly relativistic shock driven by 
the fireball into its surrounding gas. As discussed in separate chapters of this 
volume, both the spectral and temporal behavior of afterglow emission are in 
general agreement with model predictions. 

Since afterglow emission results from the interaction of the fireball with am- 
bient medium, it does not provide direct information on the evolution of the 
fireball at the earlier stage during which the GRB is produced. Nevertheless, 
afterglow observations may be used to indirectly test underlying model assump- 
tions and constrain model parameters relevant for this earlier stage. We describe 
below several afterglow observations which have important implications for the 
GRB phase of the model. 

4.1 Fireball Size and Relativistic Expansion 

Radio observations of GRB970508 allowed a direct determination of the fireball 
size and a direct confirmation of its relativistic expansion. As explained in H2.ll 
radiation seen by a distant observer originates from a conical section of the 
fireball around the source-observer line of sight, with opening angle ~ 1/-T, and 
photons which are emitted from the edge of the cone are delayed compared 
to those emitted along the line of sight, by r/2T 2 c. Thus, the apparent radius 
of the emitting cone is R — r/r{r) = 2r(r)ct where r and t are related by 
t = r /2r(r) 2 c. (A detailed calculation of fireball emission introduces only a 
small correction, R — 1.9r(r)ct |119j . Using Eq. we find that the apparent 
size of the fireball during its self-similar expansion into the surrounding medium 
is given by 

We have chosen here the normalization E = 10 52 i?52 erg since this is the energy 
inferred for GRB970508 |118I125| (note, however, the very weak dependence on 
E and n). The factor 1 + z appears due to the redshift between source and 
observer time intervals. 
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Fig. 1. The ratio between the GRB970508 afterglow radio fluxes at 4.86 GHz and 
8.46 GHz, a = log[/„(4.86GHz)//„(8.46GHz)]/ log(4.86/8.46), is shown as a function 
of time following the burst (modified from |41|1. The rapid variations at early times 
are due to narrow band diffractive scintillation, and their quenching at late times is 
due to the expansion of the source beyond a critical size given by Eq. fZQ. This is a 
direct confirmation of model predictions, according to which highly relativistic plasma 
ejection is responsible for the observed radiation. 



Scattering by irregularities in the local interstellar medium ISM may modu- 
late the observed fireball radio flux [51] . If scattering produces multiple images of 
the source, interference between the multiple images may produce a diffraction 
pattern, leading to strong variations of the flux as the observer moves through 
the pattern. In order for the diffraction patterns produced by different points on 
the source to be similar, so that the pattern is not smoothed out due to large 
source size, the apparent size of a source a redshift z = 1 must satisfy |125| 

Jt< ^,^( _^ 8 _ ) cm , m 

where v = lO^io GHz, hrs is the Hubble Constant in units of 75 km s _1 Mpc -1 , 
and the scattering measure, SM, a measure of the strength of the electron density 
fluctuations, is normalized to its characteristic Galactic value. 

Comparing Eqs. and (|24J) we find that, on time scale of weeks, the 

apparent fireball size is comparable to the maximum size for which diffractive 
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scintillation is possible. On shorter time scales, therefore, strong modulation of 
the radio flux is expected, while on longer time scales we expect diffractive scin- 
tillation to be quenched. This is exactly what had been observed for GRB970508, 
as demonstrated in Fig. ^ Observations are therefore in agreement with fireball 
model predictions: They imply a source size consistent with model prediction, 
Eq. ll'-il) which, in particular, imply expansion at a speed comparable to that of 
light. 

4.2 The Nature of the Fireball Plasma 

Due to present technical limitations of the experiments, afterglow radiation is 
observed in most cases only on time scale ^> 10 s. In one case, however, that of 
GRB990123, optical emission was detected on ~ 10 s time scale pQ. The most 
natural explanation of the observed optical radiation is synchrotron emission 
from electrons accelerated to high energy in the reverse shocks driven into fire- 
ball ejecta at the onset of interaction with the surrounding medium 81 92j, as 
explained in M3.3I This observation provides, therefore, direct constraints on the 
fireball ejecta plasma. First, it provides strong support for one of the underlying 
assumptions of the dissipative fireball scenario described in W2.2\ that the energy 
is carried from the underlying source in the form of proton kinetic energy. This 
is due to the fact that the observed radiation is well accounted for in a model 
where a shock propagates into fireball plasma composed of protons and electrons 
(rather than, e.g., a pair plasma). Second, comparison of the observed flux with 
model predictions from Eqs. I|20|) and l|21[l implies £ e ~ £b ~ 10 _1 . 

4.3 Gamma-Ray Energy and GRB Rate 

Following the determination of GRB redshifts, it is now clear that most GRB 
sources lie within the redshift range z ~ 0.5 to z ~ 2, with some bursts observed 
at z > 3. For the average GRB 7-ray fluence, 1.2 x 10 -5 erg cm -2 in the 
20 keV to 2 MeV band, this implies characteristic isotropic 7-ray energy and 
luminosity E 1 ~ 10 53 erg and L 1 ~ 10 52 erg s _1 (Here, we assume a flat 
universe with fl = 0.3, A — 0.7, and Ho = 65 km s _1 Mpc -1 . These estimates 
are consistent with more detailed analyses of the GRB luminosity function and 
redshift distribution. Mao and Mo e.g., find, for the cosmological parameters 
we use, a median GRB energy of ks 0.6 x 10 53 erg in the 50 — 300 keV band, 
corresponding to a median GRB energy of ss 2 x 10 53 erg in the 20 keV to 
2 MeV band. 

Since most observed GRB sources lie within the redshift range z ~ 0.5 — 2, 
observations essentially determine the GRB rate per unit volume at z ~ 1. The 
observed rate of 10 3 yr _1 implies i?GRB {z = 1)«3 Gpc~ 3 yr _1 (for Q = 0.3, 
A = 0.7, and Hq = 65 km s _1 Mpc -1 . The present, z = 0, rate is less well 
constrained since available data are consistent with both no evolution of GRB 
rate with redshift and with strong evolution (following, e.g., the star formation 
rate), in which Rqyib(z = 1)/Rgkb(z — 0) ~ 10 |61l68j . A detailed analysis by, 
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e.g., Schmidt |ll)lj leads, assuming -Rgrb is proportional to the star formation 
rate, to Rqrb(z = 0) ~ 0.5 Gpc~ 3 yr _1 . 

If fireballs are conical jets of solid angle Afl then, clearly, the total 7-ray 
energy is smaller by a factor AQ/Air than the isotropic energy, and the GRB 
rate is larger by the same factor. 

4.4 Fireball Geometry 

Afterglow observations suggest that at least some GRBs are conical jets, of 
opening angle 9 ~ 10 _1 corresponding to a solid angle AQ ~ 10~ 2 |26I70H 00 . 
As explained in i)2.5l the discussion in SjSJand £|3I is limited to the stage where 
the wind is ultra-relativistic, r ~ 300, prior to significant fireball deceleration 
by the surrounding medium, and is hence equally valid for both a spherical and 
a jet-like fireball. A jet like geometry has, of course, profound implications for 
the underlying progenitor (see Sj5jl: it implies that the underlying source must 
produce a collimated outflow and, if 9 ~ 10 _1 is indeed typical, it implies that 
the characteristic 7-ray energy emitted by the source is ~ 10 51 erg rather than 
w 10 53 erg implied by the assumption of isotropy. 

4.5 Fireball 7-Ray Efficiency 

Afterglow observations imply that a significant fraction of the energy initially 
carried by the fireball is converted into 7-rays, i.e., that the observed 7-ray energy 
provides a rough estimate of the total fireball energy. This has been demonstrated 
for one case, that of GRB970508, by a comparison of the total fireball energy 
derived from long term radio observations with the energy emitted in 7-rays 
|41I125| . and for a large number of bursts by a comparison of observed 7-ray 
energy with the total fireball energy estimate based on X-ray afterglow data 
|42j . Freedman and Waxman 05j demonstrated that a single measurement of 
the X-ray afterglow flux on the time scale of a day provides a robust estimate of 
the fireball energy per unit solid angle, e, averaged over a conical section of the 
fireball of opening angle 6 ~ 0.1. Applying their analysis to BcppoSAX afterglow 
data they demonstrated that the ratio of observed 7-ray to total fireball energy 
per unit solid angle, e 7 /e, is of order unity. 

The inferred high radiative efficiency implies that a significant fraction of the 
wind kinetic energy must be converted to internal energy in internal shocks, and 
that electrons must carry a significant fraction of the internal energy, i.e., that 
£ e should be close to unity. We have already shown in H3.1l and H3.2I following 
|55j . that efficient conversion of kinetic to internal energy and £ e values close 
to unity may naturally lead to ~ 1 MeV spectral break energies, in accordance 
with observations. 

5 Progenitors Clues 

The observational constraints that provide the main hints regarding the nature 
of the underlying GRB progenitors are the rapid, ~ 1 ms, 7-ray signal variability, 
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and the total 7-ray energy emitted by the source, ~ 10 53 (Z\/2/47r) erg. The rapid 
variability implies a compact object, of size smaller than a light millisecond, 
~ 10 7 cm, and mass smaller than ~ 30 M Q (the mass of a black hole with 
~ 10 7 cm Schwarzschild radius). The energy released in 7-rays corresponds to 
a 0.1(Z\/2/47r) M© rest mass energy. The most natural way for triggering the 
GRB event is therefore the accretion of a fraction of a solar mass onto a (several) 
solar mass black hole (or possibly a neutron star). The dynamical time of such 
a source, comparable to its light crossing time, is sufficiently short to account 
for the observed rapid variability and, if significant fraction of the gravitational 
binding energy release is converted to 7-rays, the source will meet the observed 
energy requirements. 

Most cosmological GRB models therefore have at their basis gravitational 
collapse of a several solar-mass progenitor to a black hole. Within the context of 
the fireball model, the observed variability is determined by the dynamical time 
of the source, which determines the variability in the ejected wind properties, 
while the GRB duration, > 10 s for long bursts, reflects the wind duration, i.e., 
the duration over which energy is extracted from the source. The characteristic 
time for the gravitational collapse is of the order of the dynamical time, i.e., 
much shorter than the wind duration. Most models therefore assume that, fol- 
lowing collapse and black hole formation, some fraction of the progenitor mass 
forms an accretion disk which powers the wind through gradual mass accretion. 
The characteristic time scale for accretion is set by the disk viscosity, which is 
uncertain and assumed to correspond to the observed GRB duration. 

Progenitor models differ in the scenario for black hole formation, and in 
the process assumed to convert disk energy into relativistic outflow. The two 
leading progenitor scenarios are, at present, collapses of massive stars |BHl28 , 
and mergers of compact objects |50I87| . In the former case, the progenitor is a 
massive rotating star, e.g., a ~ 15 M Q helium star evolved (by mass-loss) from a 
~ 30 M© main sequence star. The collapse of the progenitor's ~ 2 M© iron core 
leads to the formation of a black hole surrounded by an accretion disk composed 
of mantle plasma 1.28] . In the latter case, the merger of two neutron stars leads 
to the formation of a black hole surrounded by a disk produced by neutron star 
disruption during the merger process (see, e.g., |37j). A similar scenario involves 
neutron star disruption during a neutron star-black hole merger (see, e.g., |71p. 

Two types of processes are widely considered for the extraction of disk energy: 
neutrino emission and magneto-hydrodynamic (MHD) processes. The viscous 
dissipation of energy, driving mass accretion, heats the dense disk plasma leading 
to the emission of thermal neutrinos of all flavors. Neutrino annihilation along 
the rotation axis in the vicinity of the black hole may then produce an electron- 
positron pair plasma fireball. The fraction of rest mass which is dissipated during 
accretion is typically of order 10% (the specific energy at the last stable orbit 
of a non-rotating black- hole, at three Schwarzschild radii, is c 2 /6). Accretion 
of 0.1 M© over a second may therefore lead to a neutrino luminosity of « 
lO 52 5 erg s , of which w 1% would be deposited by annihilation to drive 
a fireball. The resulting wind luminosity, ~ 10 50 erg s _1 , may be too low to 
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drive a spherical fireball, but may be sufficient if the fireball is collimated into 
Afl/An ~ 1CT 2 . 

If equipartition magnetic fields, ~ 10 15 G, are built in the disk (e.g., by 
convective motion) the dissipated energy may be extracted electromagnetically 
from the disk. Although the process by which such a strong field is generated, 
as well as the details of the energy extraction mechanism in the presence of such 
field are not understood, there is evidence from various astrophysical systems 
(active galactic nuclei and micro-quasars, see, e.g., |74|). for the formation of 
MHD driven jets which carry ~ 10% of the disk binding energy. The presence of 
equipartition fields may also allow the extraction of energy directly, e.g., via the 
Blandford-Znajek mechanism [T^l, from the rotating black hole [72] • For rapid 
rotation, the available energy in this case is comparable to the collapsed rest 
mass. Thus, MHD processes are often invoked to drive a relativistic wind with 
efficiency much higher than that estimated for a neutrino driven wind. 

We note in this context that a possible alternative to the above models may 
be the formation from stellar collapse of a fast rotating neutron star with ultra- 
high magnetic field |1U7I1(J9| . If a fast rotating, millisecond period, neutron star is 
produced by the collapse with ~ 10 15 G field, then the resulting electromagnetic 
energy luminosity is sufficient to drive a GRB wind. 

One problem which all models are facing is the baryon loading. In order to 
allow acceleration to the high Lorentz factors implied by observations, r ~ 10 3 , 
the total mass entrained within the expanding plasma must be smaller than 
- 1CT 4 M Q , i.e., the mass-loss rate should be smaller than ~ 10 6 M Q s 1 . 
The neutrino luminosity is expected to drive mass-loss at a much higher rate. It 
is generally assumed that mass flow towards the rotation axis is inhibited (e.g., 
by high pressure of fireball plasma along the rotation axis), thus allowing the 
formation of a sufficiently baryon free fireball collimated along the rotation axis. 
In the case of a massive star progenitor, the fireball jet is assumed to form along 
the rotation axis of the star, where rapid rotation leads to lower mantle and 
envelope density. The collimation of the fireball may be due in this case to the 
presence of a low density funnel along the rotation axis, and the resulting jet 
must penetrate through the stellar mantle and envelope in order to allow the 
production of an observable GRB. Recent numerical and analytical calculations 
of the propagation of high entropy jets through stellar progenitors indicate that 
such a scenario may be viable |2I44I83| . 

Afterglow observations provide several hints which indicate that long dura- 
tion GRBs, at least those for which afterglows have been detected, are associated 
with massive star progenitors. The location of GRBs within host galaxies, the 
presences of iron lines, and the evidence for a supernova association all imply 
massive star progenitors. Since this issue is discussed in detail in other chapters 
in this volume, we address it only briefly here. 

Most GRB afterglows are localized within the optical image of a host galaxy 
|20| . This is in disagreement with simple analyses of the neutron star merger 
scenario which predict that the high velocity of such binaries should carry many 
of them outside of the host prior to merger. This result is, however, uncertain 
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since it depends on model parameters which are only poorly constrained, such 
as the distribution of initial binary separations. Evidence for the presence of 
iron lines has been found in X-ray data for two bursts |4I93| . While the presence 
of iron lines strongly suggests a massive stellar progenitor, as it indicates the 
presence of iron enriched environment, the confidence level of their detection is 
moderate. There is evidence in three cases that a supernova may be associated 
with the GRB |47I110| . The evidence is, however, not yet conclusive (see, e.g., 

124 ). Finally, the synchrotron emission produced by a shock driven by the 
fireball into its surrounding medium depends on the density of the ambient 
medium. Thus, the temporal and spectral dependence of this afterglow emission 
may distinguish between the high density environment characteristic of a massive 
stellar wind, expected to exist in the case of a massive stellar progenitor, and 
the low density ISM expected, e.g., in merger scenarios. Present observations 
are not yet conclusive, since data on time scales much shorter than one day are 
required to distinguish between the two cases [75| . 

It is clear from the above discussion, that future afterglow observations pro- 
viding more detailed information on the burst environment and location will 
play a crucial role in placing stringent constraints on progenitor models. In ad- 
dition to the points discussed in the previous paragraph, the destruction of dust 
|34I43IH6I124| and time dependence of atomic, ionic |21l7fi| and molecular H2 

34 123] lines due to photoionization may be detectable for a burst in a molecu- 
lar cloud environment characteristic of star forming regions. 

6 High Energy Protons and Neutrinos from GRB 
Fireballs 

6.1 Fermi Acceleration in GRBs 

In the fireball model, the observed radiation is produced, both during the GRB 
and the afterglow, by synchrotron emission from shock accelerated electrons. In 
the region where electrons are accelerated, protons are also expected to be shock 
accelerated. This is similar to what is thought to occur in supernovae remnant 
shocks, where synchrotron radiation of accelerated electrons is the likely source of 
non-thermal X-rays (recent Advanced Satellite for Cosmology and Astrophysics 
(ASCA) observations give evidence for acceleration of electrons in the remnant 
of SN1006 to 10 14 eV 66 ), and where shock acceleration of protons is believed 
to produce cosmic rays with energy extending to ~ 10 15 eV (see, e.g., |16| for a 
review). Thus, it is likely that protons, as well as electrons, are accelerated to 
high energy within GRB fireballs. 

We consider proton Fermi acceleration in fireball internal and reverse shocks 
(see i)2.2l and H2.4I respectively'). Since these shocks are mildly relativistic, with 
Lorentz factors .Tj — 1 ~ 1 in the wind frame (see S ^12 . 41 and 13 . 1 fl . the predicted 
energy distribution of accelerated protons is |12I16| dn p /de p cx e~ 2 , similar to 
the electron energy spectrum inferred from the observed photon spectrum. 

Two constraints must be satisfied by fireball wind parameters in order to 
allow proton acceleration to e p > 10 20 eV in internal and reverse shocks. First, 
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the proton acceleration time, t a ~ Rl/c where Rl is the proton Larmor radius, 
must be smaller than the wind expansion time |85I111I115| . td ~ fi/Tc (in the 
wind frame). This constraint sets a lower limit to the magnetic field carried by 
the wind, which may be expressed as |115j : 

W£e > 0.02r 2 2 5 e^ 20 i- 1 52 . (25) 

Here, e p = 10 20 e Pj 20 eV. Recall that £g is the fraction of the wind energy density 
which is carried by magnetic field, 4irr 2 cF 2 (B 2 / '8tt) = £,bL, and £ e is the fraction 
of wind energy carried by shock accelerated electrons. Since the electron energy 
is lost radiatively, L 1 as £ e L. 

The second constraint that should be satisfied is that the proton synchrotron 
loss time must exceed t a , setting an upper limit to the magnetic field. The latter 
constraint may be satisfied simultaneously with the lower limit to the magnetic 
field, Eq. 123), provided [HS] 

r > 130e 3 p %Arl /4 . (26) 

The constraints that must be satisfied to allow acceleration of protons to en- 
ergy > 10 20 eV are remarkably similar to those inferred from 7-ray observations: 
r > 100 is implied by observed 7-ray spectra (see S H2.ll and l2.3J) . and magnetic 
field close to equipartition, S;b ~ 1, is required in order for electron synchrotron 
emission to account for the observed radiation (see H3.1[) . 

It has recently been claimed 02] that the conditions at the external shock 
driven by the fireball into the ambient gas are not likely to allow proton acceler- 
ation to ultra-high energy. Regardless of the validity of this claim, it is irrelevant 
for the acceleration in internal shocks, the scenario considered for UHECR pro- 
duction in GRBs |111I115| . Moreover, it is not at all clear that UHECRs cannot 
be produced at the external shock, since the magnetic field may be amplified 
ahead of the shock by the streaming of high energy particles. For discussion of 
high energy proton production in the external shock and its possible implications 
see Dermer 

6.2 UHECR Flux and Spectrum 

The local (z = 0) energy production rate in 7-rays by GRBs is roughly given by 
the product of the characteristic GRB 7-ray energy, E as 2 x 10 53 erg, and the lo- 
cal GRB rate. Under the assumption that the GRB rate evolution is similar to the 
star-formation rate evolution, the local GRB rate is ~ 0.5 Gpc~ 3 yr _1 |1()1) (see 
H4.3|) . implying a local 7-ray energy generation rate of as 10 44 erg Mpc~ 3 yr -1 . 
The energy observed in 7-rays reflects the fireball energy in accelerated electrons. 
Thus, if accelerated electrons and protons carry similar energy (as indicated by 
afterglow observations 02] (see, however, |l()8j ) then the GRB production rate 
of high energy protons is 

e 2 p (dh p /de p ) z = « 10 44 erg Mpc~ 3 yr" 1 . (27) 
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Fig. 2. The UHECR flux expected in a cosmological model where high-energy protons 
are produced at a rate (<^,dh v / de p ) z= q = 0.8 x 10 44 erg Mpc -3 yr _1 as predicted in the 
GRB model (Eq. E3 solid line), compared to the High Resolution Fly's Eye (HiRes) 
|15|. Yakutsk Extensive Air Shower Array (Yakutsk) 55] and Akeno Giant Air Shower 
Array (AGASA) |105| data, la flux error bars are shown. The dashed line is the sum 
of the GRB model flux and the HiRes fit (with normalization increased by 25%) to the 
Galactic heavy nuclei component |15|. Jg oc E~ :i ' 5 , which dominates below ~ 10 19 eV. 



In Fig. El we compare the observed UHECR spectrum with that predicted 
by the GRB model. The generation rate (Eq. I27|) of high energy protons is 
remarkably similar to that required to account for the flux of > 10 19 eV cosmic- 
rays. The flux at lower energies is most likely dominated by heavy nuclei of 
Galactic origin ^j], as indicated by the flattening of the spectrum at m 10 19 eV 
and by the evidence for a change in composition at this energy (141151311461113) . 

The suppression of model flux above 10 197 eV is due to energy loss of high 
energy protons in interaction with the microwave background, i.e., to the "GZK 
cutoff' (521129) . The available data do not allow determination of the existence 
(or absence) of the "cutoff" with high confidence. The Akeno Giant Air Shower 
Array (AGASA) results show an excess (~ 2.5<7 confidence level) of events com- 
pared to model predictions above 10 20 eV. This excess is not confirmed, however, 
by other experiments. Moreover, since the 10 20 eV flux is dominated by sources 
at distances < 100 Mpc, over which the distribution of known astrophysical 
systems (e.g., galaxies and clusters of galaxies) is inhomogeneous, significant 
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deviations from model predictions presented in Fig. U for a uniform source dis- 
tribution are expected |116j . Clustering of cosmic-ray sources leads to a standard 
deviation, a, in the expected number, N, of events above 10 20 eV given by [B] 

a/N = 0.9(d /10Mpc) - 9 , (28) 

where do is the unknown scale length of the source correlation function and 
do ~ 10 Mpc for field galaxies. 

Although the rate of GRBs out to a distance of 100 Mpc from Earth, the 
maximum distance traveled by > 10 20 cV protons, is in the range of 10~ 2 to 
10 -3 yr, the number of different GRBs contributing to the flux of > 10 20 eV 
protons at any given time may be large. This is due to the dispersion, At, in 
proton arrival time, which is expected due to deflection by intergalactic magnetic 
fields and may be as large as 10 7 yr. This implies that the number of sources 
contributing to the flux at any given time may be as large as |115j ~ At x 10 -3 
yr - 10 4 . 

6.3 Neutrino Production 

A burst of ~ 10 14 eV neutrinos accompanying observed 7-rays is a natural con- 
sequence of the conventional fireball scenario |122| . The neutrinos are produced 
by 7r + created in interactions between fireball 7-rays and accelerated protons. 
The key relation is between the observed photon energy, e 7 , and the accelerated 
proton's energy, e p , at the photo-meson threshold of the Z\-resonance. In the 
observer frame, 

e 7 e p = 0.2 T 2 GeV 2 . (29) 

For r w 300 and e 7 = 1 MeV, we see that characteristic proton energies ~ 
10 16 eV are required to produce pions. The pion typically carries « 20% of the 
interacting proton energy, and this energy is roughly equally distributed between 
the leptons in the decay ir + — > /i + + — > e + + v e + + . Thus, proton 
interaction with fireball 7-rays is expected to produce ~ 10 14 eV neutrinos. 

The fraction f^(e p ) of proton energy lost to pion production is determined 
by the number density of photons in the dissipation region and is « 20% at 
high proton energy for fireball wind parameters implied by 7-ray observations 
|56I122| . Assuming that GRBs produce high energy protons at a rate given by 
Eq. (|27|l , the intensity of high energy neutrinos is |122j 

4** « 10- 9 (^) mm (l, GeV cm" 2 sr" 1 s" 1 . (30) 

Here, v x stands for v^, or v e . The neutrino flux of Eq. (|30|l is suppressed 
at high energy, > 10 16 eV, due to synchrotron energy loss of pions and muons 

gang. 

During the transition to self-similarity, high energy protons accelerated in 
the reverse shock may interact with the 10 eV to 1 keV photons radiated by the 
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accelerated electrons to produce, through pion decay, a burst of duration ~ t n of 
ultra-high energy, 10 17 - 10 19 eV, neutrinos |123| as indicated by Eq. JSJJ. The 
flux of these neutrinos depends on the density of gas surrounding the fireball. 
It is weak, and undetectable by experiments under construction, if the density 
is n ~ 1 cm -3 , a value typical for the ISM. If GRBs result, however, from 
the collapse of massive stars, then the fireball is expected to expand into a pre- 
existing wind and the transition to self-similar behavior takes place at a radius 
where the wind density is n « 10 4 enr 3 > 1 cm . In this typical 
GRB at z ~ 1 is expected to produce a neutrino fluence |30I123| 



where a = for e„ > 10 17 eV and a = 1 for e„ < 10 17 eV. The neutrino flux is 
expected to be strongly suppressed at energy > 10 19 eV, since protons are not 
expected to be accelerated to energy 3> 10 20 eV. 

6.4 Implications of Neutrino Emission 

The predicted intensity of 10 eV neutrinos produced by photo-meson interac- 
tions with observed 1 MeV photons, Eq. (|30|) - implies a detection of ~ 10 neu- 
trino induced muons per year in planned 1 km 3 Cerenkov neutrino detectors, 
correlated in time and direction with GRBs [31581122] . The predicted intensity 
of 10 17 eV neutrinos, produced by photo-meson interactions during the onset of 
fireball interaction with its surrounding medium in the case of fireball expan- 
sion into a pre-existing wind, Eq. I|31|l . implies a detection of several neutrino 
induced muons per year in a 1 km 3 detector. In this case, the predicted flux of 
10 19 eV neutrinos may also be detectable by planned large air-shower detectors 



Detection of high energy neutrinos will test the shock acceleration mechanism 
and the suggestion that GRBs are the sources of ultra-high energy protons, 
since > 10 14 eV (> 10 18 eV) neutrino production requires protons of energy 
> 10 16 eV (> 10 19 eV). The dependence of the ~ 10 17 eV neutrino flux on fireball 
environment implies that the detection of high energy neutrinos will also provide 
constraints on the GRB progenitors. Furthermore, it has recently been pointed 
out [HU that if GRBs originate from core-collapse of massive stars, then a burst 
of > 5 TeV neutrinos may be produced by photo-meson interaction while the jet 
propagates through the envelope, with TeV fluence implying 0.1 — 10 neutrino 
events per individual collapse in a 1 km 3 neutrino telescope. (The neutrino flux 
which may result from nuclear collisions in the expanding jet is more difficult to 
detect due to the low energy of neutrinos, ~ 10 GeV, produced by this process 



Detection of neutrinos from GRBs could be used to test the simultaneity of 
neutrino and photon arrival to an accuracy of ~ 1 s (~ 1 ms for short bursts), 
checking the assumption of special relativity that photons and neutrinos have 
the same limiting speed. These observations would also test the weak equivalence 
principle, according to which photons and neutrinos should suffer the same time 




(31) 
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delay as they pass through a gravitational potential. With 1 s accuracy, a burst 
at 100 Mpc would reveal a fractional difference in limiting speed of 10 -16 , and 
a fractional difference in gravitational time delay of order 10 -6 (considering the 
Galactic potential alone). Previous applications of these ideas to SN1987A, where 
simultaneity could be checked only to an accuracy of order several hours, yielded 
much weaker upper limits, of order 10~ 8 and 10 -2 for fractional differences in 
the limiting speed and time delay respectively [Hj. 

The model discussed above predicts the production of high energy muon 
and electron neutrinos. However, if the atmospheric neutrino anomaly has the 
explanation usually given [25140145] , oscillation to v T ''s with mass ~ 0.1 eV, 
then one should detect equal numbers of v^s and u T 's. Up-going t's, rather 
than /it's, would be a distinctive signature of such oscillations. Since v T 's are not 
expected to be produced in the fireball, looking for t's would be an "appearance 
experiment." To allow flavor change, the difference in squared neutrino masses, 
Am 2 , should exceed a minimum value proportional to the ratio of source distance 
and neutrino energy [Hj. A burst at 100 Mpc producing 10 14 eV neutrinos can 
test for Am 2 > 10 -16 eV 2 , 5 orders of magnitude more sensitive than solar 
neutrinos. 
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